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We investigate the branching ratios and direct CP asymmetries of i3+ — > D^vr"*" and — > D'^-k^ 
decays in the PQCD approach. All the diagrams with emission topology or annihilation topology 
are calculated strictly. A branching ratio of 10~^ and 10~^ for D^-k'^ and — > D^vr" 

decay is predicted, respectively. Because of the different weak phase and strong phase from penguin 
operator and two kinds of tree operator contributions, we predict a possible large direct CP violation: 
Af^{Bt D^-R^) « -50% and Af^{Bt -> D^vr") « 25% when 7 = 55°, which can be tested in 
the coming LHC. 

PACS numbers; 13.25.HW, 12.38. Bx 



I. INTRODUCTION 

The charmless B decays provide a good platform to 
test the Standard Model (SM) and study the CP vio- 
lation, which arouses physicists' great interest and has 
been discussed in the literature widely. But how about 
the Be decays, the h quark of which has similar property 
with that of B meson? There are some events of Be at 
Tevatron []| and will be a great number of events appear- 
ing at LHC in the foreseeable future. The progress of the 
experiments makes us to think of a question: what will be 
the theoretical prediction on the two-body non-leptonic 
Be decays? 

Different from B and Bs meson, Be meson consists of 
two heavy quarks b and c, which can decay individually. 
Because of the difference of mass, lifetime and the relative 
CKM matrix element between h and c quark, the decay 
rate of the two quarks is different, which determines the 
unique property of Be decays. Though c quark's mass 
is about one third of h quark, leading to a suppression 
of {Me/Mf,)^ , the decay of c quark can not be ignored 
because the corresponding CKM matrix element Ves is 
larger than that of b quark: Vub,Veb- Because of the 
small mass of c quark, the decay of c quark is nearly at 
non-perturbative scale, where there is great theoretical 
difficulty. Now we study the b quark decay first and leave 
the study of c quark decay to the future. 

In recent years, a great progress has been made in 
studying two-body non-leptonic B decays in perturba- 
tive QCD approach (PQCD) QCD factorization || 

and soft coUinear effective theory (SCET) 5]. Though Bj, 
decay has been studied (61] in Naive Factorization [3, Q 
many years ago, no one applies the method developed 
recently in such processes. In this paper we will use 
Be — *■ Dtt as an example to discuss the Be decays in 
the PQCD approach. 

The Be — ^ Dtt decay provides opportunities to study 
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the direct CP asymmetry. Different from the B decays 
0, Be ^ Dtt has the direct CP asymmetry even with- 
out considering the contributions from penguin opera- 
tors, because the tree contributions from the annihilation 
topology provide not only the strong phase, but also the 
different weak phase. According to the power counting 
rule of PQCD, the tree contributions from the annihila- 
tion topology is power suppressed. But the larger CKM 
matrix elements \Vcb\ enhances the contribution to make 
it larger than the penguin contributions, so the direct 
CP asymmetry of Be — > Dtt can be very large, which is 
found in our numerical analysis. 

The study of Be decay also provides opportunities to 
test kr factorization in PQCD approach. According to 
numerical analysis in the literature, the form factor con- 
tributions from Fig^usually dominate the whole decays. 
In the same way, the form factor also gives the main 
contributions in the Be — > Dtt decay according to our 
numerical analysis. Since the Be meson consists of two 
heavy quarks, the effect of kx in the Be meson can be 
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FIG. 2: Non-factorizable emission topology in Be —> Dir. 

ignored and the form factor Be ^ D only includes the fc^ 
contributions from D meson. So it is easier to study how 
important the fc^ contributions are in Be decays than 
that in the B decays because the latter need to consider 
both kx contributions of B and D meson. 

The Be — > DiT decay also provides a good platform 
to study the D meson's wave function. The D meson's 
mass Mo is not so large that it is hard to get the ideal 
wave function of D meson by the expansion of 1/Mo 
as in B meson. People use the form fitted from exper- 
imental data generally. Such discussion has been done 
by Ref. ^] in the form factor B D transition. It is 
better to discuss D meson wave function in Be — > Dtt 
for two reasons: one is that the hierarchy between Moe 
and Md {Mse ^ Md) guarantees us to apply the kr 
factorization theorem in this process, the other is that 
the wave function of Be is clean, which eliminates the 
possible uncertainty from Be meson. The experiment of 
Be decays will test how reasonable it is. As the only pa- 
rameter with large uncertainty, the wave function of D 
meson need further theoretical investigation. 




FIG. 4: Non-factorizable annihilation topology in Be Dtt. 



II. FRAMEWORK 

The hard amplitudes of these decays contain factor- 
izable diagrams (Fig. ^ , where hard gluons attach the 
valence quarks in the same meson, and non-factorizable 
diagrams (Fig.|2Jl, where hard gluons attach the valence 
quarks in different mesons. The annihilation topology 
is also included, and classified into factorizable (Fig. O 
and non-factorizable (Fig. ^ ones according to the above 
definitions. 

In the calculations of all the diagrams, we can ignore 
the kr contributions of Be meson because it consists of 
two heavy quarks. Furthermore, we can suppose the two 
quarks b and c of B^ meson to be on the mass shell 
approximately and treat the wave function of Be meson 
as S function for simplicity, so we can integrate the wave 



function Be out and the fcy factorization form turns into 
Form factor 

^ J d^ki<fD{ki)C{t)H{ki,t), (1) 
Other topology 

^ J d%d^k2<i>D{ki)^Ak2)C{t)H{ki,k2,t), (2) 

where fci(2) is momentum of light (anti) quark of D{Tr) 
meson. The non-factorizable topology includes two kinds 
of topology: emission topology (Fig. [2J) and annihilation 
topology (Fig. 21). In the above equations, we sum over 
all Dirac structure and color indices. The hard compo- 
nents consist of the hard part (H{t)) and harder dynam- 
ics (C(i)), the former H{t) can be calculated perturba- 
tively; the latter C{t) is the Wilson coefhcients which 
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runs from electro-weak scale M\y to the lower factoriza- 
tion scale t. <1>M is the wave function of D and tt meson, 
including the non-perturbative contributions in the fcy 
factorization. 

Through out the paper, we use the light-cone coordi- 
nate to describe the meson's momentum in the rest frame 
of the Be meson. According to the conservation of four- 
momentum, we get the momentum of three mesons B^, D 
and TT up to the order of r| (r2 = Md /M b^) as following: 

Mb , 
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V2 

Mb^ 

Mb^ 
V2 



(1,'^2,0t), 
(0,l-r2,OT), 



(3) 



where we have neglected the small mass of pion and 
higher order terms of r2- Such approximation will be 
used in the whole paper. 



III. CALCULATION OF AMPLITUDES 



A. Wave Function 

Be meson consists of two heavy quarks such that the 
small Aqcd can be ignored (Aqcd = Mbc — Mb — Ale <C 
Mc or Mb), so as the quark transverse momentum kx- 
In principal there are two Lorentz structures in the B 
or Be meson wave function. One should consider both 
of them in calculations. However, it can be ar gue d that 
one of the contribution is numerically small [lOj, thus 
its contribution can be neglected. Therefore, we only 
consider the contribution of one Lorentz structure, such 
that we can reduce the number of input parameters 



iB^+MBj75S{x^Me/MBj. 

The other two mesons' wave functions read: 



l5{>P2+MD)Mx,b), 

75 f'acjy^ix) + Mo^j5K{x) 




(4) 



(5) 



(6) 

where iVc = 3 is color degree of freedom, and AIqtt — 
M^/{mu + md), n_ = (0,1, Ot) oc P3, n+ = (1,0, Ot), 

g0123 _ 2 

The momentum fraction of the light quark in the 
three mesons can be defined as: xi = ke/ Pb^, X2 = 
^2 I P2 = / ■ tti6 Be meson, there are also an- 
other relation between xi and = Mb/MB^- xi+rb = 1. 



B. Effective Hamiltonian 

The effective Hamiltonian for the flavor-changing b 
d transition is given by [T]| 



H. 



oft 



Gf 
V2 



i=3 



(7) 



with the Cabibbo-Kobayashi-Maskawa (CKM) matrix el- 
ements Vq — VqdV*f^ and the operators 

o['^^ = {d^qj)v-A{qjbi)v-A , 



O^'^-' = {d^qi)v-A{qjbj)v-A 



)v~A 



)V-A 



03 = {dibi)v-A^{qjqj)\ 

q 

04 = {dibj)v-A^{qjqi)\ 

q 

05 = {dibi)v-A^{qjqj)v+A 

q 

{dib.j)v-A^{qjqi)v+A 
q 

l^idzbiW-A^^egiqjqj) 



Os = ^{d,bj)v-A^eq{qjqi) 

3 _ 
O9 = T^i^zb 



V+A , 



V+A 



^i)V~A 



V-A 



3 _ ^ 

Olo = -{d,bj)v-A^eq{qjq 



i}V-A 



(8) 



with i and j being the color indices. Using the unitary 
condition, the CKM matrix elements for the penguin op- 
erators O3-O10 can also be expressed as Vu + Ve — — Vt- 
The Be Dtt decay rates have the expressions 



327r ' ' 



(9) 



The decay amplitude A oi Be^ Dtt process from all the 
diagrams can be expressed in the following: 



D°TT + 



VuU.Fl + Mji) + VeifBFj + Mj) 
+fBFr + fBFr + Mr+Mr' 



Af- (10) 



V2A 



= VuiUF^f^ + MJ2) - Ve{fBFj + Mj) 



'Vt {LFP' + UF^^ + LF^' 
+M^2'+Mf,' + M^,^-fBF^ 

-Ibf; 



PI 



P3 



PI 



MP) 



(11) 



where F{M.) denotes factorizable (non-factorizable) am- 
plitudes, the subscript e(a) denotes the emission (annihi- 
lation) diagrams. The subscript 1(2) denotes the process 
B+ D°7r+ {B+ D+Tr°), the superscript T{P) de- 
notes amplitudes from the tree (penguin) operators, and 
/sc (/tt) is the Be (tt) meson decay constant. The de- 
tailed expressions of these amplitudes are shown in Ap- 
pendix^ 

From ea. H10lll|) . we can see that unlike i?^, i?°(i?°) 
decays, we have three kinds of decay amplitudes with 
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different weak and strong phases: penguin contributions 
proportional to Vt and two kinds of tree contributions 
proportional to Vc and 14, respectively. The interference 
between them gives large direct CP violation which will 
be shown later. 

As stated in the introduction, the two diagrams in 
Fig.l give the contribution for Be ^ D transition form 
factor, which is defined as 



We calculate in PQCD and get: 
4/s 



(12) 



'TV, 



nl nOO 

■.ttCfMI^ / dx2 / b2db2 
" Jo Jo 



4>d{x2M) 



X I [2n -X2- in - 2x2)r2 + {x2 - 'i.nyl] 

+ [{l-xi)r2{2-r2)] 

X as {tl^^)SD (4'^ )He2 (ae ,l3e2,b2)}ll3) 

which is the similar expression as F^^ without Wilson co- 
efficients in the appendix. The numerical results of i^-f 
can be found in Tabled the form factor F+ is 0.169+°;"^ 
including the uncertainty of (£n, which is comparable 



with previous calculations 



OlUJn, 



F+ 



0.40GeV 0.45GeV 0.50GeV 



0.154 



0.169 



0.174 



TABLE I: Form factor _F+ in the different values of dd 



C. Input parameters 

For D meson wave function, two types of D meson 
wave function are usually used in the past literature: one 
is i 



4>d{x) 



X exp 



fDxil-x){l + aDil-2x)} 
1 



-(LUDbY 



(14) 



in which the last term 
kr distribution; the other 



' e xp [-^( ^ 
aer HI 111 i 



LUobY], represents the 



(koix) 



/2N, 



:fDx{l~x){l + aD{l-2x)}, (15) 



which is fitted from the measured B Dtv decay spec- 
trum at large recoil. The absence of the last term in the 
Ea. (|14|l is due to the insufficiency of the experimental 
data. 

Though the wave function of D meson turns more com- 
plicated when it runs at a velocity of about 0.6c, the light 
quark's momentum must be less than because the 

mass of c quark is by far larger than Aqcd : A/c 3> Aqcd , 
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FIG. 5: D meson wave functions: the dashed line for case 1 
&2, the solid line for case 3. 



so the wave function of D meson should be strongly sup- 
pressed in the region X2 > 1/2 even the kx contribu- 
tions are considered. In order to satisfy the above con- 
dition, we give up the D wave functions above and con- 
struct a new wave function, which also fits the measured 
B Dlv decay spectrum at large recoil: 



{x,b) 



Nd [x{l 
1 
2 



exp 



x)f X 
xMd 



-^{-ofb' 



where Nu is a normalization constant to let 



4>D(x,b) = — 1^= 



(16) 



(17) 



The behavior of all D meson wave function can be seen 
in the FigEl Our choice of the third case has a broad 
peak at the small x side, which characterize the mass 
difference of nic and m^. 

The TT wave functions 0, we adopt are calculated 
by QCD sum rules and shown in the Appendix IbI 

The other input parameters are listed below 11^. IT^: 

fs, = 480 MeV, fo = 240 MeV, U = 131 MeV, 
WD=0.45GeV, Mo^ = 1.60GeV, = 0.3, 
Mb, = 6.4 GeV, Mb = 4.8 GeV, 
Md = 1.869 GeV, = 170 GeV, 



Mw = 80.4 GeV, tb± = 0.46 x lO^^^s, 
Gf = 1.16639 X 10"^ GeV"^. 

The CKM parameters used in the paper are: 

Vub 



0.085 ± 0.020, 

Vcb 

\Vcb\ = 0.039 ± 0.002, 



R 



\_}?I2_ 
A 



ub 



(18) 

(19) 
(20) 
(21) 



The CKM angle 03 = 7 is left as a free parameter to 
discuss CP violation, defined by 



arg 



arg(K\)- 



(22) 
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D. Numerical analysis 

We fix 7 = 55° to discuss the central value of numerical 
results first. 

Both process — > D°7r+ and _D+7r° are tree- 

dominated. The branching ratios and main contributions 
are give in Table IHI from which we can see that the 
branching ratio of i?+ D^tt^ is much larger than that 
of D^Ti^ . Though they are both tree-dominated 

process, their branching ratios and percentage of differ- 
ent topologies in the whole process are obviously differ- 
ent. Because the annihilation topology give the same 
contributions to both processes, despite a \/2 factor, the 
difference only comes from the emission topology. In 
the process 5+ — > Z?°7r+, contributions from factorizable 
emission topology dominate the whole tree contributions 
for the large Wilson Coefficients C2 -I- Ci /iVc in Ea. (|JT3)l . 
which occupy about 93% of total even when the effect of 
CKM matrix element is considered. (|A„| < |Ac|). On the 
contrary, contributions from factorizable emission topol- 
ogy in the process D~^n^ are suppressed because 
the Wilson Coefficients Ci and C2/NC in Ea. (|A14|l cancel 
each other approximately. From the Table Hll we also find 
that contributions from factorizable annihilation topol- 
ogy are at the same order of non-factorizable emission 
topology. 



23.0 
-0.379 + 0.863i 
-3.35 + 5.49i 
2.52 - 1.92j 



0.763 
0.854 - 2.16i 
-3.35 + 5.49i 
2.52 - 1.92i 





P 














10% 


40% 












Br 


0.978 X 10"" 


0.196 X 10"*^ 



TABLE II: Brancfi ratios and main contributions from tree 
operators (lO'^GeV ). 

The ratio of the penguin contributions over the tree 
contributions is about 10% in the process — > Z3°7r+ 
and about 40% in the process B+ D+tt° (Table HHl. 
The reason for the difference is the following: the term 
2r^^P(a;3) in the Ea. (|A4|l from Oq,Os operator contri- 
butions, having no factors like 2:3 to suppress its integral 
value in the end-point region and leading to large en- 
hancement compared with other penguin contributions. 
But the most important reason is that the tree contri- 
bution is suppressed in the process B^ £)+7r° due to 
small Wilson coefficients Ci -f C2/3 but not suppressed 
in the process B+ D^TT+. The Oq,Os. contributions 
also affect the dependence behavior of the branching ra- 
tio and the direct CP asymmetry on the CKM angle 7 in 
the process Bf^ D^n^ , which will be discussed in the 
following. 

The correlation between Br{B^ Dtt) and r^r is 
shown in Fig^l Because twist— 3 terms of tt wave func- 
tion do not contribute to the form factor fEa. (jAllA2p ). 
the variation of r^r affects the process i?+ D^tt^ more 
heavily than the process i?+ — + D'^tt^ , where the latter 
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FIG. 6: The correlation between Br{Bc Dtt) and r,r 



dominated by the Be ^ D form factor diagrams. When 
= 1.4, the twist-3 contributions are about 25% in the 
process i?+ D^tt"*". In the process B^ D^n^, the 
twist-3 contributions with a relative minus sign cancel 
some of the twist-2 contributions. When r-^ = 1.4, the 
branching ratio of 5+ _D+7r° is about four times of the 
branching ratio with only twist-2 contributions. When 
= 0, the twist-3 contributions vanish and only the 
contributions from twist-2 terms in the tt wave function 
are left. The corresponding branching ratio is reduced to 
0.95 X 10-6 in the process B+ D"n+ and 0.092 x 10"^ 
in the process B^ D^tt^ respectively. 





B+ ^ D\+ 


B+ D+7T° 


LOD 


= 0.40GeV 


1.03 


0.128 


UlD 


= 0.45GeV 


0.978 


0.196 


UJD 


= 0.50GeV 


1.19 


0.199 



TABLE III: Branch ratios in the unit 10 for different ujd ■ 

As the only free parameter with large uncertainty, the 
value of lud is the key point to the whole prediction in 
the calculations of Be Dtt. In the Table ITTll we discuss 
the branching ratio in three groups of different to values: 
LVD = 0.40GeV, UJD = 0.45GeV and ivd = 0.50GeV, from 
which we see that the variation of ujo affects the process 
B+ D^TT+ slightly, but affect the process B+ D+tt° 
heavily. 

According to the CKM parametrization shown in the 
Ea. p9l22|l . the decay amphtudes of Be Dtt can be 
written as: 



VuTu + VeTe ~ VtP 

1 Te 



P) 



VuiTu + P) 1 



1 - 



P 



RTa + P 



(23) 
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where z 
(5 — ars 



1 












R 










T^ + P 



and the strong phase 



f Tc + P \ 
\T^+P J ' 



J, I , from our PQCD calculation the numer- 
ical value of which is 0.28 and 123° for B+ D°tt+, re- 
spectively. The emission topology in this channel is only 
about one time larger than the annihilation topology due 
to the smah CKM factor \Vu/Vc\. 
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FIG. 7: The correlation between the averaged branching ratio 
and 7 in the process D'^tt*. 



The corresponding conjugate decay of i3+ — s- Dtt 
reads: 



v:{Tu+p) 



1 — ze 



and the averaged branching ratio for B- 
reads: 



(24) 



Br = \{\M\^ + \Mt) 

- i|K(T„ + P)|2 [l-2zcos7COS(5 + z2] ,(25) 

which is the function of CKM angle 7. Its numerical 
result depends on 7 significantly: the larger 7, the smaller 
the averaged branching ratio, because cos 5 < 0. The 
explicit correlation between the averaged branching ratio 
Bf D°{D^)tt^ and 7 is shown in the Fig0 
The direct CP violation A'^" is defined as 

dir^ |Af(B+^g°( + ). + (°)|^-|M(B-^g°'-).-(")|^ 
|M(B+^D0(+)7r + (0)|" + |Af(B-^D0(-)7r-(")|" ' 

There are two different tree contributions and one kind 
of penguin contribution with different strong and weak 
phases, which will contribute to the CP asymmetry. Us- 

idir 



ing ea. (|23l24|) . A^}^ can be simplified as 



A 



dir 



2zsini5sin7 



1 — 2z cos 5 cos 7 



(27) 



which is proportional to sin 7 approximately. This is 
shown in Fig|Hl When 7 = 55°, the direct CP asym- 
metry is about —50% in the process i?^ Z?°7r+. 

The Be Z?+7r" becomes a little more complicated: 
the tree contributions from the emission topology M]^ 
(in Table I) is suppressed due to the small Wilson co- 
efficients Ci -I- C2/3. In this case, the three different 
contributions with different weak and strong phases (two 








zo 40 60 80 100 120 140 1 eoyi 80 

FIG. 8: The correlation between the direct CP violation and 
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FIG. 9: The correlation between the averaged branching ratio 
and 7 in the process B^ D^n°. 



tree contributions and one penguin contributions) are at 
the same order of magnitude. We can still use Ea. H25|) 
and Ea. H27|l to get the behavior of the branching ratio 
and the direct CP asymmetry on 7. Now the numerical 
values of z and 6 are 3.1 and —20° respectively. Dif- 
ferent from the averaged branching ratio of the process 
Be D+tt", the averaged branching ratio of the process 
Be — > D^7T~^ becomes smaller when 7 becomes larger be- 
cause cos (5 > 0. The behavior of the branching ratio and 
the direct CP asymmetry does not change much, but the 
shape of the former turns more sharp. In one word, the 
branching ratio of Bf D^tt" shown in Fig^are more 
sensitive to the value of 7, which is quite different from 
the case for B^ D^n^, but the direct CP asymmetry 
of Bf — > D*7r° shown in Fig|Sldoes not change greatly 
because the large uncertainty from 7 cancels in the ratio 
of the direct CP asymmetry. When 7 = 55°, the direct 
CP asymmetry is about 25% in the process B^ D^tt'^. 
As pointed out in ref.^^, the CP asymmetry is sensitive 
to the next-to-leading order contribution, which is more 
complicated, the result shown here should be taken care- 
fully. 



IV. CONCLUSION 

In this paper we discuss the process Be D°TT+ and 
Be D'^tt'^ in the PQCD approach and get their branch- 



ing ratios imtoM X 10" 



6 and 1.96^^;^^ x 10" 



respec- 
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tively. We also predict the possible large direct CP asym- work is partly supported by National Science Foundation 
metry in the two processes: A'^p{B^ D^n^) « —50% of China, 
and (^^ ^ D^TT°) « 25% when 7 = 55°. The pos- 
sible theoretical uncertainty are also analyzed. We hope 
it can be tested in the coming experiments at Tevatron, 
LHC and the super-B factory. 
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APPENDIX A: CONTRIBUTIONS FROM ALL THE DIAGRAMS 

1. Contributions from factorizable diagrams 

All diagrams are sorted into two kinds: emission topology and annihilation topology shown in Fig lll3l and Fig l2l4l 
The factorizable tree contributions from emission topology read: 

= -^I^TtCfMI^ f dx2 f h2dh2<t>D{x2M) 

V 2jVc Jo Jo 

X { [2n -X2- in - 1x2)r2 + {X2 - 2rb)r2] El^^^^'"^\t^^'^)H,^{a,,(i,M 

+ [(1 - x^)r2(2 ~ r2)\ El^^'''^\tf>)H,2{a,,fi,2. 62)}, (Al) 

f2 = -^^rK^CpMl^ [ dx2 [ h2dh2C^D{x2M) 

V^J^^c Jo Jo 

X { [-2 + + (1 - An + X2)r2 + [n - 2x2 + 2)rl] £^{1^^^ )H,i {ae,Pei , ^2) 

- [xi + 2(1 - 2xi)r2 - (2 - a;i)r2] £;f^3(^(2))^^^(^^^ /3e2, 62)}. (A2) 

Because b and c are both massive quarks, there is no coUinear divergence in the Be D transition, so the threshold 
resummation needn't to be considered. In all the expressions, T denotes the contributions from tree operators, PI 
denotes the penguin contributions with the Dirac structure {V — A)®{V — A)^ P2 denotes the penguin contributions 
with the Dirac structure {V — A) ^ {V + A), and P3 denotes the penguin contributions with the Dirac structure 
[S — P) (S + P); the subscript e(a) denotes the factorizable emission (annihilation) diagrams, the subscript ne [no) 
denotes the nonfactorizable emission (annihilation) diagrams. 

The factorizable tree contributions from annihilation topology read: 

Fj(^l)=8^Cj^M|^ I dx2dx3 f b2db2b3db3cj)D{x2,b2) 

Jo Jo 

X { [{x3 - (1 + 2x3)r2)0,(x3) + r2r,((l + 2x3)0^(2:3) - (1 - 2x3)^(2^3))] 

xE]^^P^\t(^^)Ha{aa,Pal,b2,b3)St{x3) 

- [x2(l-r2)<^,(x3)+2r2r,(l+X2)C(x3)] E^^P^\t^^^)Ha{aa, /3a2,b3,b2)St{x2)} , (A3) 

FP^~IQttCfMI^ I dx2dx3 / b2db2b3db3(t)D{x2,b2) 
" Jo Jo 

X { [-r2<?!>^(x3) + - X3 + (2 + X3)rl)4,l + r^X3(l - rl)(j)l{x3)\ 

- [x2r24>Ax3) + 2r^{l^{l^X2)rl)^l{x3)\ EP{t^;^'i)Ha{aa, Pa2,b3,b2)St{x2)], (A4) 
where the factor St{x) is the jet function from the threshold resummation [20| : 

St{x)= ^r(/^J) \ x{l-x)r, c = 0.3. (A5) 
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The factors eJ^^^ (t) contain the Wilson coefficients a(t) at scale t and the evolution from t to the factorization scale 
1/6 in the Sudakov factors S{t): 



E. 



where (t) , S'tt (t) , the Sudakov factors, are defined as 



Soit) 



s,{t)^s{x^p-,b^) + s{{i~x^)p^M) + 'i / 



(A6) 

(A7) 
(A8) 



and s{Q,b) is given as |21 



i/b M 



^(27.-l-ln2) + CWn^^"^(^) 



67 TT^ 
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nf H — /?n In 

9 3 27 3 2 



(A9) 



where the Euler constant je = 0.57722 • • •, and — — as/vr is the quark anomalous dimension. 
The hard functions H are 



Hei{a,P,b) 
He2{a,P,b) 



Kojab) ~ Ko{(3b) 
13^ -a^ 



-Koiab), 



(l-a;i)(a;i - rj)' 



(AlO) 

(All) 
(A12) 



where Ko,Io,Ho and Jo are the Bessel functions of order 0. It is implied that the transformed Bessel functions Kq 
and Iq become the corresponding Bessel functions with real variable when their variables are complex. 
The Wilson Coefficients read: 



,P3 



Ci + — + Clo + — , 



(A13) 



N, 



C5 



C4 



C3 



Nr. 



(A14) 



C9 



C7 



Nr. 



a2it) = -{C, + ^ 



Cio 



Nr 



Clo 



Nr 



(A15) 



All the Wilson coefficients Ci above should be evaluated at the appropriate scale t. The hard scale t's are chosen as 
the maximum of the virtuality of internal momentum transition in the hard amplitudes, including 1/6^: 

= max(|ae|,|/3ei|,l/62), 

= max(|ae|,|/3e2|,l/62), 

= max (1/3,11,1/62,1/63), 

= max (1/3,21,1/62, 1/^3), 
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where 



/3f 



ol 



o2 



(1-xi - .X2)(.xi -r|)M: 



Bel 



(1 



-2:2(1- 
Xi)(xi 



= -X2X3.Ml^{l - rl 
= ~x^Ml^{l-rl), 
= -X2Ml^{l-rl). 



(A16) 



2. Contributions from non-factorizable diagrams 

Different from factorizable diagrams, non-factorizable diagrams include convolution of all three wave functions and, 
of course, the convolution of Sudakov factors. Their amplitudes are: 



M, 



T(P1) 



M. 



M. 



X 
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-'^c Jo Jo 

x| [1 - - cca - (1 - - a;2)r2 - {x2 - '^.x^yl] El'f^\t^^'')Ha{ane,(inelMM) 

+ [(2a;i + a;2 - 0:3 - 1) + (1 - a;i - a;2)r2 + (-2a;i -X2+ 2xi)rl] El^/^\t^;^^) /3ne2, 62, 63)}, 

(A17) 

—Trr^CFfsMl / (ia;2da;3 / h2db2hidhi(j)D{x2,h2) 
Nc Jo Jo 

X { [(1 - xi - 0:3 + (2 - 2x1 - X2 - x:i}r2 + (1 - xi - X2 + x^)rl)(l)P^{x:i) 

+ (1 - xi - X3 + {x2 - X'i)r2 + (-1 + xi^X2+ x^)rl)(t>l{xi)\ E^"^ ^{t^^^) Ha{ane, Iinel,b2,bz) 
+ [{xi - X3 + (2x1 + X2 -X3- l)r2 + (xi + X2 + .X3 - 2)rl)(f)P^{x3) 

+ ( - a;i + X3 + {X2 + X3- l)r2 + {xi + X2 - X3)rl)<pl{x3)] ^^^e Ha{a„e,l3ne2, 62, 63)}, (A18) 

/■I POO 

TtCf/bMI^ / dX2dX3 / h2dh2h3dh34)D{X2,h2)4>'K{X3) 

Jo Jo 

I [-2 + 2x1 + X2 + a;3 + (1 - xi - a;2)r2 + (2 - 2xi -X2- 2x3)rl] £;^e^(^i^^) i?a(ane, /3„ei, &2, ^s) 

+ [-Xl +X3 + {Xl +X2- l)r2 - {X2 + 2X3 " 2)r^] E^^ ^{t'^^^) Ha{ane, f3ne2, ^2, 63)}, (Al9) 

-TTT^CFfEM^ I dX2dX3 / 62rf&2&3ci&3?!'r>(a;2, ^2) 
Nc Jo Jo 

x{[( -Xl + X2 + r2)^7r(a;3) + (-2x1 + a;2 + X3 + 4r2)r2r^<^(x3) + (x2 - a;3)r2r^C(x3)] 

X E^el^^ {tl^l ) Hna {ana , l^nalM) 

+ [(1 - Xl - X3 - Tfe + (-X2 + 2x3 + ''fe)rf)^^(x3) + (2 - 2x1 - a;2 - X3 - 4r6)r2r^?!>^(x3) 

Hx2-X3)r2r^4.l{X3)]El^^'\t(^l)Hna{ana,Pna2M)}, (A20) 

-^■kCfIbM^^ / dx2dx3 / h2dh2b3dh3(t)D{x2,b2) 

^^c Jo Jo 

X I [{-Xl +X2- r2)r2 (t>T,ix3) + [xi - X3 + r2 + (xi - X2 + X3)r|)r„^(x3) 

+ (xi +X2 - X3 + (-Xl +X2 +X3)r2)r^C(X3)]£;^g^l(4a^)if„a(a„a,/J„ai,62) 

- [(-1 - rb + Xl + X2)r2(^„(x3) + (1 + rfc - Xl - X3 + (1 + rfc - Xl - X2 + X3)r2)r^ ^^^(a^s) 

+ (1 + rb - Xl - X3 + (-1 - r6 + Xl + X2 + X3)rl)rMx3)]Kel{t^nl) Hna{ana,Pna2,b2)}. (A21) 

where the hard kernel H^a is defined as 

Ko{ab) - Kom 



Hna{a,l3,b) = 



(A22) 



10 



and the factor E{t) turns into: 



(A23) 



where the Wilson coefficients a read: 



,P2 



it) = 
(t) = 
(t) = 



C3 + , 

C2, 



P2 



P3 



it) - -C5 



2^10, 



2^8 



(A24) 



The hard scale t's are chosen as the maximum of the virtuality of internal momentum transition in the hard amplitudes, 
including l/bf. 

= max(|a„e|, |/3nel|, I/&2, 1/&3) , 

tf' 

= max (|a„e|, |/3ne2|, 1/^2, 1/^3) , 

t^^^ ^ max(|Q;„a|, |/3„ai|, 1/62) , 

t^^^ = max (|Q;„a|, |/3„a2|, 1/62) , 



where 



al = {\ - xi - X2)ixi - rl)Ml^, 

Pnei - -il-xi-X2)[il-X3){l-ri)--xi+rl]Ml^, 

Pne2 = -il-xi-X2)[x3{l-rl)-xi+rl]Ml^, 

al = -X2X3Ml^{l-rl), 

final = Xi [X2 + X^ll - rl)] , 

/3«ai = (l-xi) [x2+ 0:3(1 -r^)] Af|^. 



(A25) 



APPENDIX B: THE tt MESON WAVE FUNCTIONS 

The different distribution amplitudes of tt meson wave functions are given as |l5l Il6j | 

(l)^{x) = —=f„x{l-x) l + 0.44Cf^(22;-l) + 0.25Cf^(2a;-l) 
V6 L 



€i^) 
Kix) 



u 

2V6 



1 + 0A3Cl^^(2x - 1) + 0.09C]^^(22; - 1) 



(1 - 2x) [1 + 0.55(10a;2 - lOx + 1)] 



with the Gegenbauer polynomials 



Cl^^it) - liSt^ - 1), C'j^t) = i(35t4 - 30f + 3), 



C2^^it) = ^(5i" - 1), C^'"it) = ^(2K^ - ue + 1) 



3/2/ 
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(Bl) 
(B2) 
(B3) 



(B4) 
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